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ABSTRACT 

In this paper, we consider a multihop wireless sensor network 
(WSN) with multiple relay nodes for each hop where the amplify- 
and-forward (AF) scheme is employed. We present a strategy to 
jointly design the linear receiver and the power allocation parame- 
ters via an alternating optimization approach that maximizes the sum 
rate of the WSN. We derive constrained maximum sum-rate (MSR) 
expressions along with an algorithm to compute the linear receiver 
and the power allocation parameters with the optimal complex am- 
plification coefficients for each relay node. Computer simulations 
show good performance of our proposed methods in terms of sum 
rate compared to the method with equal power allocation. 

Index Terms — Maximum sum-rate (MSR), power allocation, 
multihop, wireless sensor networks (WSNs) 

1. INTRODUCTION 

Recently, there has been a growing research interest in wireless sen- 
sor networks (WSNs) as their unique features allow a wide range of 
applications in the areas of defence, environment, health and home 
|T| . They are usually composed of a large number of densely de- 
ployed sensing devices which can transmit their data to the desired 
user through multihop relays |2j. Low complexity and high energy 
efficiency are the most important design characteristics of commu- 
nication protocols (3i and physical layer techniques employed for 
WSNs. The performance and capacity of WSNs can be significantly 
enhanced through exploitation of spatial diversity with cooperation 
between the nodes |2|. In a cooperative WSN, nodes relay signals to 
each other in order to propagate redundant copies of the same signals 
to the destination nodes. Among the existing relaying schemes, the 
amplify-and-forward (AF) and the decode-and-forward (DF) are the 
most popular approaches ||4l . 

Due to limitations in sensor node power, computational capac- 
ity and memory |T), some power allocation methods have been pro- 
posed for WSNs to obtain the best possible SNR or best possible 
quality of service (QoS) (5] at the destinations. The majority of the 
previous literature considers a source and destination pair, with one 
or more randomly placed relay nodes. These relay nodes are usually 
placed with uniform distribution J6], equal distance 0, or in line (8) 
with the source and destination. The reason of these simple consider- 
ations is that they can simplify complex problems and obtain closed- 
form solutions. A single relay AF system using mean channel gain 
channel state information (CSI) is analyzed in [9j, where the outage 
probability is the criterion used for optimization. For DF systems, a 
near-optimal power allocation strategy called the Fixed-Sum-Power 



with Equal-Ratio (FSP-ER) scheme based on partial CSI has been 
developed in (6). This near-optimal scheme allocates one half of the 
total power to the source node and splits the remaining half equally 
among selected relay nodes. A node is selected for relay if its mean 
channel gain to the destination is above a threshold. Simulations 
show that this scheme significantly outperforms existing power al- 
location schemes. One is the 'Constant-Power scheme' where all 
nodes serve as relay nodes and all nodes including the source node 
and relay nodes transmit with the same power. The other is the 'Best- 
Select scheme' where only the node with the largest mean channel 
gain to the destination is chosen as the relay node. 

In this paper, we consider a general multihop wireless sensor 
networks where the AF relaying scheme is employed. Our strategy 
is to jointly design the linear maximum sum-rate (MSR) receiver 
(w) and the power allocation parameter (a) that contains the optimal 
complex amplification coefficients for each relay node via an alter- 
nating optimization approach. It can be considered as a constrained 
optimization problem where the objective function is the sum-rate 
(SR) and the constraint is a bound on the power levels among the 
relay nodes. Then the constrained MSR solutions for the linear re- 
ceiver and the power allocation parameter can be derived. The pro- 
posed strategy and algorithm are not only applicable to simple 2-hop 
WSNs but also to general multihop WSNs with multi relay nodes 
and destination nodes. Another novelty is that we make use of the 
Generalized Rayleigh Quotient 1131 to solve the optimization prob- 
lem in an alternating fashion. 

This paper is organized as follows. Section 2 describes the mul- 
tihop WSN system model. Section 3 develops the joint MSR re- 
ceiver design and power allocation strategy. Section 4 presents the 
proposed alternating optimization algorithm to maximize the sum 
rate. Section 5 presents and discusses the simulation results, while 
Section 6 provides some concluding remarks. 

2. SYSTEM MODEL 

Consider a general m-hop wireless sensor network (WSN) with mul- 
tiple parallel relay nodes for each hop, as shown in Fig. 1 . The WSN 
consists of No source nodes, N m destination nodes and N r relay 
nodes which are separated into m — 1 groups: N\,N2, ... ,-/V m _i. 
We will focus on a time division scheme with perfect synchroniza- 
tion, for which all signals are transmitted and received in separate 
time slots. The sources first broadcast the No x 1 signal vector s to 
the first group of relay nodes. We consider an amplify-and-forward 
(AF) cooperation protocol. Each group of relay nodes receives the 
signal, amplifies and rebroadcasts them to the next group of relay 
nodes (or the destination nodes). In practice, we need to consider 



the constraints on the transmission policy. For example, each trans- 
mitting node would transmit during only one phase. In our WSN 
system, we assume that each group of relay nodes transmits the sig- 
nal to the nearest group of relay nodes (or destination nodes) directly. 



Cooperative 
relay nodes 
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nodes 




Fig. 1. m-hop WSN with Nq sources, N m destinations and N r relays. 



Let H a denote the N\ x No channel matrix between the source 
nodes and the first group of relay nodes, denote the N m x N m -\ 
channel matrix between the (m — l)th group of relay nodes and 
destination nodes, and Hj-i,, denote the Ni x N-i channel matrix 
between two groups of relay nodes as described by 
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Ni is a row 



ij — [h s ,j,i, h s j,2, —, h s ,j,N ] for j = 1,2 
vector between source nodes and the jth relay of the first group of 
relay nodes, h m -i,j = [h m -i,j,i,h m -i,j,2,...,h m -i t j t N m _- L ] 
for j = 1, 2, ...,N m is a row vector between the (m — l)th 
group of relay nodes and the jth destination node and hi_i,j = 
[hi-ij t i, hi-i,j,2, —, fti-ij^-i] for j = 1, 2, Ni is a row 
vector between the (i — l)th group of relay nodes and the jth relay 
of the ith group of relay nodes. The received signal at the ith group 
of relay nodes (Xi) for each phase can be expressed as: 
Phase 1: 



Phase 2: 



Xi = H s s + Vi, 
Yi = Fixi, 

X2 = Hi, 2 Aiy 1 + v 2 , 
y 2 = F 2 x 2 , 



(2) 
(3) 

(4) 
(5) 



Phase i: (i = 2, 3, m — 1) 

Xi = H i _i,iAi_iy^_ 1 + 

At the destination nodes, the received signal can be expressed as 
d = H d A m _iy m _ 1 + Y d , (8) 



(6) 
(7) 



where v is a zero-mean circularly symmetric complex additive white 
Gaussian noise (AWGN) vector with covariance matrix <j 2 I. A; = 
diagjai.i, 6^2, «i,JVi} is a diagonal matrix whose elements rep- 
resent the amplification coefficient of each relay of the zth group. Fj 
denotes the normalization matrix which can normalize the power of 
the received signal for each relay of the ith group of relays. (see the 
appendix to find the expression of Fj.) Please note that the property 
of the matrix vector multiplication Ay = Ya will be used in the next 
section, where Y is the diagonal matrix form of the vector y and a 
is the vector form of the diagonal matrix A. At the receiver, a linear 
MMSE detector is considered where the optimal filter and optimal 
amplification coefficients are calculated. The optimal amplification 
coefficients are transmitted to the relays through the feedback chan- 
nel. And the block marked with a Q[-] represents a decision device. 

3. PROPOSED JOINT MAXIMUM SUM-RATE DESIGN OF 
THE RECEIVER AND THE POWER ALLOCATION 

By substituting (2)-(7) into (8), we can get 

d =Cn, m _lS + Cl, m _lVl + C2,m-lV2 
+ ■•■ + C m _i m-lVm-1 + v d 



— Co,m-lS + y Ci, m -lVi + Vd 



where 



and 



Ci.j — 



I, if i > j. 



(9) 
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B, 



Bo = H s (11) 

Hi,i + iAiF« for i = l,2, m - 2 (12) 

B m ~i = H £ iA m _iy m _ 1 (13) 

We focus on the system which consists of one source node. There- 
fore, the expression of the Sum Rate (SR) for our m-hop WSNs is 
expressed as 
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■log 2 



1 + 
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where w is the linear receiver, and ( ) H denotes the complex- 
conjugate (Hermitian) transpose. Let 



and 



4> — Co,m-lC(). m _l 



Z — ^ ^ Ci,m — iCi - 



(15) 
(16) 



The expression for the sum-rate can be written as 
1 , /„ a'i w H chy/\ 1 



SR = - log 2 1 + 
m 

where 
and 



log 2 (l + ax) (17) 



(18) 



(19) 



w H Zw 

Since — log 2 (l + ax) is a monotonically increasing function of x 
(a > 0), the problem of maximizing the sum rate is equivalent to 



maximizing x. In this section, we consider the case where the total 
power of the relay nodes in each group is limited to some value Pr.i 
(local constraint). The proposed method can be considered as the 
following optimization problem: 



[v/opt, ai, opt, ...,a m -L„ p t] = arg max 



w H (frw 



',»1,..,>™-1 w H Zw' (20) 

subject to Pi = Pt.i, i = 1, 2, m — 1. 

where Pi as defined above is the transmitted power of the ith group 
of relays, and Pi = Ni+iafai. We can notice that the expression 
"jj^" in (20) is the Generalized Rayleigh Quotient, therefore the 
optimal solution of our maximization problem can be solved: w pt is 
any eigenvector corresponding to the dominant eigenvalue of Z~ <p. 

In order to obtain the optimal power allocation vector a op t, we 
rewrite w and the expression is given by 



Table 1. Summary of the Proposed Algorithm 
Initialize the algorithm by setting 

Ai = \JW^ 1 for 1 = l > 2 ' 771 - 1 
For each iteration: 

1. Compute </> and Z in (15) and (16). 

2. Using QR algorithm or power method to compute the 
dominate eigenvector of Z _1 </>, denoted as w op t. 

3. For i = 1, 2, m — 1 

a) Compute M t and N t in (24) and (25). 

b) Using the QR algorithm or power method to compute the 
dominate eigenvector of Mj, denoted as a^. 

c) To ensure the local power constraint a.f opt a.i, op t = N T ' Z , 



compute a i>opt 



N i + l s 'i a i 



w g </>w afM,a ; . 

— 77 = —77 ; — 77 ; 77 , lor 1 = 1,2,..., m—i, 

w H Zw af diagjwf PjWijai + wf TjWj 

(21) 

where 

Mi = diag{wf Ci+i, m _iHi 1 i+iFi}Co,<_iCo[ i _ 1 x 
diag{lfHf 1+1 Cg. lim _ 1 w<}, 

Pi = Cj-|_i jm _iHi l j-|_iFi}(2^^,_ 1 Cfc l j_iC^_ 1 )x 

diag{Fi Hii+jCi+x^^!} 

and 

Ti = (Z~2k = i + 1 Cfe,m-lC fe . m _l). 

Since the multiplication of any constant value and a eigenvector 
is still the eigenvector of the matrix, we can express the receive filter 
as 



opt 



(22) 



opt 



Therefore, we can obtain 



m AT H 
Mi ( 2^ Cfc, m _iC fc>m _i)Wi = 1 = (23) 



Pt., 



By substituting (23) into (21) and using Mi and Nj given above to 
simplify the expression of (21), we obtain 



w H <£w afM 4 ai 



w H Zw afNiai 



fori = 1, 2,..., m- 1 (24) 



We can notice that the expression * H N y* in (26) is the Generalized 

Rayleigh Quotient, therefore the optimal solution of our maximiza- 
tion problem can be solved: Sn, pt is any eigenvector corresponding 
to the dominant eigenvalue of Mi, and satisfying &^opt&i,opt = 
ff-^ ■ The solutions of w op t and ai, op t depend on each other. There- 
fore it is necessary to iterate them with an initial value of a; (i = 
1, 2, 77i — 1) to obtain the optimum solutions. 



4. PROPOSED ALTERNATING MAXIMIZATION 
ALGORITHM 

In this section, we devise our proposed alternating maximization al- 
gorithm which computes the linear receive filter and the power allo- 
cation parameters that maximize the sum rate of the WSN. In partic- 
ular, we employ two methods to calculate the dominant eigenvectors. 



The first one is the QR algorithm 1181 which calculates all the eigen- 
values and eigenvectors of a matrix. We can choose the dominant 
eigenvector among them. The second one is the power method H 1 31 
which only calculates the dominant eigenvector of a matrix. There- 
fore, the computational complexity can be reduced. Table 1 shows 
a summary of our proposed algorithm which will be used for the 
simulations. 

Tong, we should include a short paragraph here discussing the 
required complexity and the convergence issues. 

5. SIMULATIONS 

In this section, we numerically study the sum-rate performance of 
our proposed joint MSR design of the receiver and power allocation 
methods and compare them with the equal power allocation method 
[ 6 1 which allocates the same power level equally for all links from 
the relay nodes. We consider a 3-hop (m=3) wireless sensor net- 
work. The number of source nodes (No), two groups of relay nodes 
(Ni, N2) and destination nodes (-/V3) are 1, 4, 4, 2 respectively. We 
consider an AF cooperation protocol. The quasi-static fading chan- 
nel (block fading channel) is considered in our simulations whose el- 
ements are Rayleigh random variables (with zero mean and unit vari- 
ance) and assumed to be invariant during the transmission of each 
packet. In our simulations, the channel is assumed to be known at the 
destination nodes. For channel estimation algorithms for WSNs and 
other low-complexity parameter estimation algorithms, one can refer 
to 1151 and 1161 . During each phase, the sources transmit the QPSK 
modulated packets with 1500 symbols. The noise at the destination 
nodes is modeled as circularly symmetric complex Gaussian random 
variables with zero mean. When perfect (error free) feedback chan- 
nel between destination nodes and relay nodes is assumed to trans- 
mit the amplification coefficients, it can be seen from Fig. 3 that our 
proposed method can achieve better sum-rate performance than the 
equal power allocation method. When using the power method to 
calculate the dominant eigenvector, it can get a very similar result to 
the QR algorithm. In practice, the feedback channel can not be error 
free. In order to study the impact of feedback channel errors on the 
performance, we employ the binary symmetric channel (BSC) as the 
model for the feedback channel and quantize each complex ampli- 
fication coefficient to an 8-bit binary value (4 bits for the real part, 
4 bits for the imaginary part). Vector quantization methods 1171 can 
also be employed for increased spectral efficiency. The error prob- 
ability (Pe) of BSC is fixed at 10 -3 . The dashed curves in Fig. 3 
show the performance degradation compared with the performance 



when using a perfect feedback channel. To show the performance 
tendency of the BSC for other values of Pe, we fix the SNR at 10 dB 
and choose Pe ranging form to 1CP 2 . The performance curves are 
shown in Fig. 4, which illustrates the sum-rate performance versus 
Pe of our two proposed methods. It can be seen that along with the 
increase in Pe, their performance becomes worse. 



where 

f a 2 s \h a ,j\ 2 +al, fori = l, 

y for i = 2, 3, m. 

(26) 
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Fig. 2. Sum-rate performance versus SNR of our proposed joint maximum 
sum-rate design of the receiver and power allocation strategy for a 3-hop 
WSN, compared with equal power allocation method. 
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- Local Constraint, Powor Method (BSC 8Bits) 
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Fig. 3. Sum-rate performance versus Pe of our proposed joint strategy when 
employing BSC as the model for the feedback channel. 



6. CONCLUSIONS 

A joint MSR receiver design and power allocation strategy has been 
proposed for general multihop WSNs. It has been shown that our 
proposed strategy achieves a significantly better performance than 
the equal power allocation method. Possible extensions to this work 
may include the study of the complexity and the requirement for the 
feedback channel. 



7. APPENDIX 

Here, we derive the expression of F; which is denoted in Section 2. 



Fi(cr?H s Hf + o£l)Ff , for i = 1, 

F i pi_ ll< A i _ 1 .B(y i _ 1 y£ 1 )A£. iHf_ M + <x£l]Ff 

for i = 2,3, m. 

(27) 
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